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ABSTRACT 

Observed X-ray spectra of hot gas in clusters, groups, and individual galaxies are commonly fit with a single- 
temperature thermal plasma model even though the beam may contain emission from components with different 
temperatures. Recently, Mazzotta et al. pointed out that thus derived r spec can be significantly different from 
commonly used definitions of average temperature, such as emission- or emission measure-weighted T, and 
found an analytic expression for predicting r spec for a mixture of plasma spectra with relatively hot temperatures 
(T > 3 keV). In this Paper, we propose an algorithm which can accurately predict r spec in a much wider range 
of temperatures (T > 0.5 keV), and for essentially arbitrary abundance of heavy elements. This algorithm can 
be applied in the deprojection analysis of objects with the temperature and metallicity gradients, for correction 
of the PSF effects, for consistent comparison of numerical simulations of galaxy clusters and groups with the 
X-ray observations, and for estimating how emission from undetected components can bias the global X-ray 
spectral analysis. 

Subject headings: X-rays: galaxies: clusters — X-rays: galaxies — methods: N-body simulations — tech- 
niques: spectroscopic 



1. INTRODUCTION 

Temperature of the hot gas filling the volume of galaxy 
clusters and groups is the primary diagnostic of properties and 
physical processes in these objects. An incomplete list of ap- 
plications includes the study of radiative cooling and feedback 
mechanisms in the cluster centers; distribution of heavy ele- 
ments in the intracluster (ICM) and interstellar (ISM) media; 
estimation of the cluster mass either through the virial M — T 
relation or application of the hydrostatic equilibrium equation 
(e.g.. lEvrard et alJl996t[Sarazinll988HMathewsll978l) . 

ICM temperature is usually measured by fitting its observed 
X-ray spectrum. Generally, the spectrum is integrated within 
a beam which contains several components with different T 
and metallicity. Current detectors, such as CCDs onboard 
Chandra and XMM-Newton, cannot spectrally separate emis- 
sion from different components. Also, statistical quality in the 
vast majority of cases is insufficient to detect the presence of 
several emission components in the total spectrum. Therefore, 
single-temperature models are commonly fit to the integrated 
spectrum with the hope that the derived T is a representative 
average of the temperatures within the beam. 

Recently, Mazzotta et al.] J2004I) pointed out that the tem- 
perature derived from the X-ray spectral analysis, r spec , is 
significantly different from commonly used averages, such as 
the emission measure-weighted Tem (volume-averaged with 
weight w = p 2 ), or emission-weighted Te {w = p 2 A(T), 
where A(T) is the plasma emissivity per unit emission mea- 
sure). The "spectroscopic" temperature is biased towards 
lower-temperature components and it is generally lower than 
either of Te or Tem- 

As discussed below, r spec depends not only on the input 
spectrum but also on the energy dependenc e of the ef fective 
area of the X-ray telescope. However, Mazzotta et al. (2004) 
were able to find a simple analytic weighting scheme which 
predicts T spec for a known distribution of plasma tempera- 
tures and is sufficiently accurate for both Chandra and XMM- 
Newton, as long as the minimum temperature is suff iciently 
high, r m ; n > 3 keV. This work can be applied (Mazzotta et al. 



2004) for realistic comparison of the cluster numerical simu- 
lations and observations, for estimating detectability of hydro- 
dynamic phenomena (shocks) in the X-ray data, for modeling 
the 3D distributions in the presence of temperature gradients, 
and for estimating how the presence of undetectable compo- 
nents can bias t he global cluster properties inferred from the 
X-ray analysis (Rasia et al. 2005). 

Unf ortunately, the weighting scheme of iMazzotta et alJ 
(2004) can be applied only for relatively high temperatures, 
Tmin k. 3 keV, because it was developed for continuum- 
dominated spectra. It is important to be able to accurately 
predict r spec in the lower- temperature regime. For example, 
many of the cool clumps within the clusters, whose pres- 
ence biases the global r spec , probably have temperatures typ- 
ical of galaxy groups or i ndividual galaxies, T ~ 1 - 2 keV 
( iMotl et : al JI20041 iNagai et ai]l2003t iDolag et al.ll2004 . An- 
other application is in the analysis of the cluster regions out- 
side half the virial radius where the ICM temperature drops 
below 50% of its peak value near the center (see Vikhlinin et 
al. 2005 for recent results on the cluster temperature profiles). 
An algorithm to predict r spec for low-temperature plasmas is 
required for the X-ray analysis of low-mass clusters, the ob- 
jects which will provide the bulk of cosmological constraining 
power in the forthco ming Sunyaev-Zeldovich effect surveys 
(lHaiman et all2001l) . 

In this Paper, we develop an algorithm which accurately 
predicts r spec in the entire range of X-ray temperatures (T > 
0.5 keV), and for nearly arbitrary range of plasma metallici- 
ties. The algorithm is successful because it explicitly accounts 
for the low-energy line emission as the primary temperature 
diagnostics for low-T plasmas. The general plan is as fol- 
lows. In §|3 we briefly review how temperature is derived 
from the X-ray spectral analysis. Extreme cases of purely 
line-dominated and continuum-dominated spectra are consid- 
ered in §|3] Spectroscopic temperature for in the case of real- 
istic metallicities, Z = 0.1 — 2 Solar, can be predicted using 
the results for these extreme cases, as discussed in §|4] 
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Fig. 1. — Top: Observed Chandra spectrum of a T = 0.8 keV galaxy group 
USGC S152. Histogram shows the MEKAL model spectrum convolved with 
the detector response. Bottom: Model spectrum binned into 10 eV energy 
channels. 

2. X-RAY SPECTRAL ANALYSIS 

Before proceeding to discussion of the temperature- 
averaging techniques, we briefly review some technical as- 
pects of the ICM temperature determination by current X- 
ray telescopes. X-ray spectrum of plasma enriched by heavy 
elements is superposition of continuum bremsstrahlung and 
emission lines. Advanced codes exist which can com- 
pute spectra of collisional-dominated plasma in a broad 
range of temperature s and metallic ites: the Raymond-Smith 
model dRavmond & Smith! Il977h. MEKAL dMewe et al.1 
119851 i Kaastra & Mewell993tlLiedahl et al.ll995l) . and APEC 
( Smit h et al.l2001l) . We use the MEKAL model in this Paper; 
the results for other codes should be nearly identical. 

Observed spectrum is significantly distorted because of the 
energy-dependent effective area of the X-ray telescopes and 
finite detector spectral resolution. For rigorous mathe matical 
model of the instrumental response to input spectra, see lDavisl 
J2001I) . X-ray CCDs of the ACIS camera onboard Chandra 
and EPIC onboard XMM-Newton have a ~ 100 eV energy res- 
olution (FWHM). Individual emission lines are blended when 
observed with such resolution and line and continuum emis- 
sion at low energies cannot be fully separated (Fig^. 

Since it is impossible to reconstruct the input spectrum di- 
rectly, its parameters are determined from fitting a model to 
the observed spectrum. Typically, the free parameters are the 
overall normalization, temperature, and abundance of heavy 
elements. The most temperature-sensitive features in the ob- 
served spectrum are the location of high-energy exponential 
cutoff in the continuum, if the temperature is sufficiently low 
so that the cutoff is within the instrument bandpass; the over- 
all slope of the continuum emission; the location of the low- 
energy emission line complex (see below). The best-fit model 
tends to describe the most statistically significant of these fea- 
tures. Which feature takes precedence is determined by both 
the input spectrum and instrument characteristics such as rela- 
tive sensitivity at low and high energies, the instrument band- 
pass, and the background. Therefore, different instruments 
will, generally, measure different T if the input spectrum con- 
sists of multiple temperature components. 
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Fig. 2. — Average energy for emission lines in the MEKAL spectra con- 
volved with the Chandra BI response, as a function of plasma temperature. 
The lines from bottom to top are for N H = 0, 4 X 10 20 , 10 21 , and 2 X 10 21 . 

In this paper, we focus on the single-temperature fits with 
free abundance of heavy elements, performed with the Chan- 
dra back-illuminated (BI) CCDs in the 0.7-10 keV energy 
band. We then show that the results are nearly identical in 
the same energy band for the Chandra front-illuminated (FI) 
CCDs, and only slightly different for the XMM-Newton and 
ASCA detectors (§0. More serious modifications will be re- 
quired for very different instrumental setup, e.g. for high- 
resolution X-ray spectrometers or very different energy bands. 

3. r spec IN EXTREME CASES 
3.1. Average Temperature from Fitting Emission Lines 

Ionization states of heavy elements and ion excitation rates 
sensitively depend on the plasma temperature. This makes the 
relative stre ngths of emission lines a n attractive temperature 
diagnostics (Sarazin & Bahcall 1977). The brightest lines in 
the typical ICM and ISM spectra are those of O VII and O VIII 
at E ~ 0.6 keV for T < 0.5 keV, the iron L-shell complex at 
E ~ 1 keV for T = 0.5 - 3 keV, and Fe XV and Fe XVI 
lines at 6.7 and 6.95 keV for T > 3 keV. In practice, the 
high energy Fe lines are not important for temperature deter- 
mination because at the relevant temperatures, the spectrum 
is continuum-dominated. The Fe L-shell complex and O VII 
and O VIII lines, on the contrary, are very important because 
the plasma emission at low temperatures is line-dominated for 
commonly found metallicities (Z > 0.1 Solar). 

The individual iron L-shell lines are not resolved by X-ray 
CCD detectors and the complex is observed as a single broad 
bump in the spectrum (Fig^. As the plasma temperature in- 
creases, the complex is shifted towards higher energies. In 
fact, its location is the strongest temperature diagnostics in 
the low-T regime. This suggests the following proposition: 
when a single-temperature model is fit to the line-dominated 
spectrum, the best-fit T is such that the average energy in the 
input and model spectra are equal. The average energy can be 
defined as 

(E) = f(T) = (1) 
H^- 
where Si is the observed count rate in channel i and £, is 
the nominal energy corresponding to this channel. The count 
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Fig. 3. — Single-temperature fits to line-dominated spectra. The input spectra consist of two components with temperatures r m ; n and 7" max , and with relative 
emission measures / m j n and (1 - / m i n ). Left: results for (T m i R , T max ) = (0.25, 0.5), (0.5, 1), and (1,2) keV. Right: results for mixtures of 0.5 and 1.5 keV, and 1 
and 3 keV components. On both panels, points show T for a single-temperature fit to the simulated spectra. Solid line shows predictions of eq. lTld . Dotted lines 
correspond to naive weighting, <T> = / min r min + (1 - / min ) T max . 



rates s, depend on the temperature, interstellar absorption, and 
detector sensitivity as a function of energy. Figure |3 shows 
(E) as a function of temperature assuming that the observa- 
tion is performed with the Chandra BI CCDs. 

Obviously, the average energy for a mixture of components 
with total count rates S j and temperatures Tj is given by 



<£>tot = 



(2) 



iMazzotta et alJ (12004 . Mazzotta et al. noted that the spec- 
trum of the high-temperature bremsstrahlung can be approxi- 
mated by a linear law, s(E) « a(T)—b(T)E, within the limited 
bandpass of Chandra and XMM-Newton and this leads to the 
following weighting scheme for computing the average tem- 
perature, 

J v wTdV 



(T) = 



f v wdV ' 



The proposition formulated above implies that a single- 
temperature fit to the combined spectrum can be predicted by 
inverse function of eq. fT). 



where 



w = p 2 T~ a , 



a a 0.75. 



(4) 



(5) 



fit, lines 



(3) 



The functions (E) = f{T) can be pre-computed and tabulated 
and then eq. Q provides an easy and fast method for predict- 
ing T for a single-temperature fit. 

The a ccuracy of th is technique can be tested by direct 
XSPEC (lArnaudl 19961) simulations (Fig.E). The single-F fit 
can be accurately predicted for a mixture of 2 components 
with different temperatures as long as the dynamical range is 
not too large, T m;ix /T m i n < 2. For larger temperature differ- 
ences, predictions become less accurate because the emission 
complexes are well-separated and the composite spectrum is 
bimodal. The best fit in this case tends to model the brighter 
component and ignore the weaker one. Still, eq. provide 
qualitatively correct predictions for the single-temperature 
fit which are much more accurate than weighting by emis- 
sion measure (right panel in Fig|3J- In more realistic cases, 
the temperature is distributed continuously between T m i D and 
Tmax- The composite spectrum then will be unimodal and pre- 
dictions for 7fit should be quite accurate even when Tm^/T^n 
is large. 

3.2. Average Temperature for Purely Continuum Spectra 

Let us now consider the opposite case of spectra with 
zero metallicity. Predictions for the single-temperature fit to 
the continuum-dominated spectra were recently derived by 



Mazzotta et al. demonstrated that this formula is remarkably 
accurate for both Chandra and XMM-Newton, as long as all 
temperatures are sufficiently high, T > 3.5 keV. 

We would like to extend the weighting scheme |3] into 
the lower temperature regime. The obvious problem here 
is that the weights |5) are strongly skewed towards lower- 
temperature components. For T — » 0, w — > oo but in reality 
it should become zero because the spectrum is shifted below 
the bandpass of the X-ray detectors. A heuristic modification 
of eq. could then be 



c(T)p 2 T~ a , 



(6) 



where c(T) is the detector sensitivity to bremsstrahlung spec- 
tra characterized by the observed photon count rate within the 
energy band of interested for a spectrum with unit emission 
measure. For T — > 0, c(T) — > and thus suppresses the 
weights for the low- temperature components. 

Surprisingly, we find that eq. (6) with a = 0.875 works 
accurately in a very broad temperature range. The results 
of XSPEC simulations for mixtures of two components with 
7max/7min = 4 and T m i n = 0.25, 0.5, 1, 2, and 4 keV are 
shown in Fig@] The approximations by eq. 1416 1 are shown 
by the solid lines. They are accurate in these cases to within 
AT < 0.15 keV or AT/T < 4%, whichever smaller. Approx- 
imations using eq. J5] (dotted lines) are equ ally accurate at 
high te mperatures but fail for T < 2 keV, as IMazzotta et alJ 
(2004) warned. 
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Fig. 4. — Single-temperature fits to continuum-dominated spectra consist- 
ing of two components with temperatures r m j n and 7" max , and with relative 
emission measures f m [ n and (1 - / m i n ). Simulations were performed for 
Chandra BI CCDs and for (T min ,T mlIX ) = (0.25,1), (0.5,2), (1,4), (2,8), 
and (4, 16) keV. Filled circles show best-fit T for a single-temperature mode l 
fit in the 0.7-10 keV energy band. Approximations for r spcc using 14161 . 
Approximations using eq. |5] are shown by dotted lines for comparison. 

The drawback of our modification of the Mazzotta et al. 
weighting scheme is that it is no longer purely analytic. The 
function c(T) should be pre-computed and tabulated because 
it is unique for each observation setup, which includes the 
telescope effective area as function of energy, Galactic ab- 
sorption, and energy range used for spectral fits. Since the 
setup can be quite arbitrary, c(T) should be implemented as a 
computer code. However, the gains in accuracy and tempera- 
ture range warrant these complications. 

4. T spec FOR SPECTRA WITH TYPICAL 
METALLICITIES 

Thermal emission of ICM and ISM with typical metal- 
licities, Z = 0.1-1 Solar, is usually not purely line- or 
continuum-dominated. To predict the single-temperature fit 
for such metallicities, we need to find a valid combination of 
results for the extreme cases considered above. 

A possible approach is suggested by results of XSPEC sim- 
ulations shown in Fig. [5] The filled circles in this Figure cor- 
respond to single-temperature fits to mixtures of T = 1 and 
3 keV spectra with metallicities of Z = 0.1 Solar, and of 
T — 2 and 6 keV with Z = 0.3 Solar. The approximations 
for the line-dominated (Z —> oo) and continuum-dominated 
(Z = 0) regimes derived in §|3]are shown by dotted and dashed 
lines, respectively. For both cases, T\i ne < T cont and the real 
single-temperature fit, r spec , is between these two regimes. 
Also note that r spec approaches 7"ii ne for large values of fm m , 
when the contribution of the line emission to the total count 
rate becomes more significant. The total spectrum is more 
continuum-dominated for small values of / m i n , and we ob- 
serve that r sp ec approaches r cont . This suggests that the pa- 
rameter x, 

^spec ^line 
^cont — ^line 

should depend on the relative contribution of the line and con- 
tinuum emission to the total flux in the energy band used for 
spectral fitting. 



Fig. 5. — Single-temperature fits to spectra with realistic metallicities. The 
input spectra have temperatures of 1 and 3 keV and metallicity Z = 0. 1 Solar 
(lower curves), and T = 2 and 6 keV and Z = 0.3 Solar (upper curves). Filled 
circles show best-fit T for a single-temperature model fit in the 0.7-10 keV 
energy band. Approximations for continuum- and line-dominated spectra 
(see §|5J are shown by dashed and dotted lines, respectively. 




Fig. 6. — The values of parameter x in eq. as a function /ij ne , fraction of 
the line emission in the total flux in the 0.7-10 keV energy band. Data points 
were obtained from XSPEC simulations spanning a wide range of tempera- 
tures, r max /r m i n = 2, 3, and 4, r m j n = 0.5, 1, 2, and 3, and metallicities of 
0.1, 0.3, and 1 Solar. Solid line shows analytic approximation of eq.l8l. 

We have performed a large number of XSPEC simulations 
of two-component spectra for a wide range of r m ; n , r max , and 
metallicites. Figure |6] shows the values of parameter x as a 
function of /ii, le , fraction of the line emission in the total flux 
in the 0.7-10 keV energy band which we use for spectral fit- 
ting. The data points in this Figure were derived from simula- 
tions with Tmax/rmin = 2, 3, and 4, r m ; n = 0.5, 1, 2, and 3, and 
metallicities of 0.1, 0.3, and 1 Solar. Even though these cases 
probe very different regimes, all x seem to follow a nearly 
universal function, which can be approximated as 
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Fig. 7. — Single-temperature fits to multiple-component spectra with T 
distributed between r m i n and T mdx and emission measure Vp 2 oc T y (see 
text). Filled and empty circles show the simulation results for T m i n = 1 keV 
and T max = 6 keV, and for T m m - 0.5 keV and T mllx = 3 ke V, respec- 
tively. Approximations using the algorithm summarized in § 14. H are shown 
by solid lines. Temperatures weighted by emission measure and using eq. |5] 
are shown by dotted and dashed lines, respectively. Blue points show the 
temperatures measured by the XMM-Newton EPIC-pn detector for the same 
input spectra, and blue solid line are the predictions of our model for this 
detector. 

with p = 1, Ai = 0.19, and A 2 = 0.25 (solid line). As ex- 
pected, x — > 1 and r spec — > T mnt for /i; ne — » 0, and x — > and 
r sp ec — * T\i ne for /ii ne — > 1 . The transition between the two 
regimes occurs for /jj ne 0.2. Apparently, emission lines at 
this point become the strongest feature in the observed spec- 
trum and therefore they are the main driver for the single- 
temperature fit. 

To further test the universality of eq. JH) we performed sim- 
ulations for different values of Galactic absorption, Nh = 0, 
10 20 , 5xl0 20 and2xl0 21 cm 2 . The values of x obtained from 
these simulations were within the scatter of the data points in 
Fig-El This shows that at least for the same instrument setup 
(effective area and energy band of the spectral fit), the pre- 
scription for combining continuum- and lines-based temper- 
atures is universal. Therefore, to properly combine r cont and 
Time, we need to know only the fraction of the line emission in 
the total flux, fu n e(T), a quantity which is easy to pre-compute 
and tabulate. 

4.1. Putting It All Together: The Algorithm 

To efficiently implement the method outline above, the fol- 
lowing functions should be pre-computed and tabulated: 

c(T) — observed photon count rate for purely continuum 
spectra with unit emission measure, as a function of tempera- 
ture; 

A(T) — observed photon count rate for line emission for 
spectra with unit emission measure and Solar metallicity; 

(E) = f(T) — average energy of the line emission (eq. (T|). 

The continuum-based tempera ture for the composite spec- 
trum, r cont , is obtained using eq. 14 161 and the total continuum 
flux is given by the integral 

^cont = j c(T) P 2 dV. (9) 

The total flux and mean energy of the line emission are given 



TABLE 1 

Parameter values for Chandra, XMM-Newton, and ASCA detectors 



Parameter 


Chandra 


XMM-Newton 


ASCA 


BI, FI 


PN 


MOS 


SIS GIS 


a (eq. El) 


0.875 


0.790 


0.900 


0.875 0.790 


/SCeq-II) 


1.00 


0.75 


1.00 


0.80 0.75 


Ai (eqTgl) 


0.19 


0.27 


0.19 


0.20 0.26 


A 2 (eq. |U) 


0.25 


0.225 


0.22 


0.22 0.30 



by 

^line = j A(T)Zp 2 dV, (10) 

(E) = F^ e ff(r ) A(T)Zp 2 dV, (11) 

where Z(x, y, z) is the metallicity distribution. The line-based 
temperature, T\i 1K , is given by eq. Q. To properly combine 
Tcont an< ^ ^ line ' we mst com P ute me fraction of the line emis- 
sion in the total observed flux, fn ne = i* , ii ne /(Fii ne + F cont ), 
then find the value of x from eq. and finally compute the 
spectroscopic-like temperature as 

Tspec — x r cont + (1 — X) T\i ne . (12) 

To test how this scheme works in realistic cases, we per- 
formed XSPEC simulations for five-component spectra with 
temperatures equally log-spaced between r m ; n and r max , 

T t = r mm (r max /r mm ) !/4 , i = o, 1,2,3,4, (13) 

and emission measures distributed as a power law of temper- 
ature, 

V iP 2 KTj. (14) 

The relative contribution of the low- and high-temperature 
components to the overall spectrum is controlled by the value 
of index y, which we vary in the range from y = -1 to 
y — 3. Results of XSPEC simulations for such complex 
spectra are shown in Fig.0 The results of single-temperature 
fits in the case of r m i n = 0.5 keV and r max = 3 keV, and 
7min = 1 keV and r max = 6 keV are shown by open and filled 
circles, respectively. Solid lines show the predicted values of 
r spec . Clearly, we are able to predict the single-temperature 
fit very accurately. The residuals in Fig.^are less than 7% 
or 0.05 keV, whichever smaller. A similar high accuracy of 
predictions is found in all other realistic cases we checked. 
Our algorithm becomes inaccurate only in extreme cases — 
for example, when the input spectrum has two components 
with similar flux and very large temperature difference. Such 
cases are easily identifiable in practice because the single- 
temperature model provides a very poor fit to the data. 

To illustrate the gain in accuracy achieved by our algo- 
rithm, we also computed weighted temperatures for the spec- 
tra shown in Fig.[7] using weighting by emission measure 
(w = p 2 ) and by w given in eq. 0- The results for these cases 
are shown in Fig.0by dotted and dashed lines, respectively. 

5. RESULTS FOR XMM-NEWTON AND ASCA 

In the discussion above, we have used XSPEC simulations 
performed for the Chandra BI CCDs. We now should check 
how sensitive are parameter values in eq. ]6) and JS] to the 
choice of the X-ray detector. The complete analysis was re- 
peated for Chandra FI CCDs, and also for the XMM-Newton 
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and ASCA detetors, using the 0.7-10 keV energy band for 
spectral fitting in all cases. We find no significant difference 
in the results for the Cha ndra BI and FI CCDs. However, the 
parameters in equations 16181 derived for XMM-Newton and 
ASCA are slightly different. This is expected because these 
instruments have a different relative effective area at low and 
high energies, and hence difference sensitivities to thme con- 
tinuum emission. The parameters for all detectors are listed 
Tabled Blue lines in Fig.Qshow predictions of our model 
for the XMM-Newton observations. 



projection analysis of objects with the temperature and metal- 
licity gradients, consistent comparison of numerical simula- 
tions of galaxy clusters and groups with the X-ray observa- 
tions, and estimating how emission from undetected compo- 
nents can bias the global X-ray spectral analysis. 

The algorithm requires precomputed tables of several pa- 
rameters of the observed spectra as a function of temperature. 
Fortran code which implements these computations is publi- 
cally available from the following WEB page: 
http : //hea-www . harvard . edu/~ alexey/mixT. 



6. CONCLUSIONS 

We presented an algorithm for predicting results of single- 
temperature fit to the X-ray emission from multi-component 
plasma. The algorithm is accurate in a wide range of temper- 
atures and metallicities. Possible applications include the de- 
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